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Calculation of nuclear spin-dependent parity-nonconserving amplitude for 

(7s, F = 4) -> (7s, F = 5) transition in Fr 
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Many-body calculation of nuclear spin-dependent parity-nonconserving amplitude for 7s, F = 
4 — » 7s, F = 5 transition between hyperfine sublevels of the ground state of 211 Fr is carried out. 
The final result is (7s, F = 5||dpNc||7s, F — 4) = —0.49 x 10 -10 i K a.u., where k is the dimensionless 
coupling constant. This is approximately an order of magnitude larger than similar amplitude in 
Cs. The dominant contribution to k is associated with the anapole moment of the nucleus. 

PACS numbers: 32.80Ys, 11.30.Er, 31.30.Jv 



a. Introduction. In this work we calculated nuclear 
spin-dependent parity-nonconserving (PNC) amplitude 
for 7s, F = 4 — ► 7s, F = 5 transition between hyper- 
fine structure components of the ground state of the odd 
isotope of francium 211 Fr. Three effects contribute to this 
amplitude m : the interaction of an electron with the nu- 
clear anapole moment (AM), the electron- nuclear neutral 
current interaction, and the combined action of the nu- 
clear spin-independent electron-nucleus weak interaction 
and the hyperfine interaction. 

The AM a was introduced by Zel'dovich g just after 
the discovery of parity violation. A first realistic model 
for the AM of the nucleus was suggested in Refs. || Q. 
There it was shown that for heavy nuclei a ~ A 2 ' 3 , where 
A is the number of nucleons. AMs of the nuclei with un- 
paired proton are expected to be few times larger than 
for the case of unpaired neutron. Because of that for 
atoms with large and odd Z the AM contribution to 
the spin-dependent part of the PNC amplitudes domi- 
nates over that of the electron- nucleon neutral currents. 
The third contribution is also ~ A 2 / 3 , but is numerically 
smaller g, [| (see Eq. (g) below). Note that the neu- 
tral current and hyperfine contributions to the nuclear 
spin-dependent PNC amplitude are well known from the 
standard model. Therefore, any measurement of the re- 
spective coupling constant k will give unambiguous in- 
formation about AM of the nucleus. 

For the optical transitions in heavy atoms the spin- 
independent PNC amplitudes are approximately two or- 
ders of magnitude larger than the spin-dependent ones. 
Because of that the AM was measured experimentally 
only for cesium M. This measurement provided a valu- 
able probe of the relatively poorly understood parity non- 
conservation in nuclei [|3| |9|. Further experimental and 
theoretical investigations of AM are very important both 
for nuclear physics and for physics of the fundamental 
interactions. 

An alternative possibility to observe the spin- 
dependent PNC amplitudes was suggested in jiqj : in 
the rf transitions between the hyperfine components of 



the ground state of an atom the spin-independent am- 
plitudes are negligible and the dominant PNC effect is 
caused by the AM. Using the rf resonator one can have 
an additional enhancement of the PNC effect by placing 
the gas cell in the node of the magnetic and the antin- 
ode of the electric rf fields @. The PNC effect can be 
also enhanced in the strong dc magnetic field ||12[ . The 
new cooling and trapping techniques allow to increase 
the intensity of the rf transitions making these experi- 
mental schemes much more realistic. At present there is 
an ongoing project of measuring PNC effects in francium 
p~3|, [hi and the observation of the AM in the hyperfine 
transition can be a valuable addition to this project. 

Semi-empirical calculations of nuclear spin-dependent 
amplitudes for transitions between hyperfine sublevels of 
the ground state were already carried out for Cs and Tl 
Jiof and for K and Cs [jll). Fr is the heaviest of alkali- 
metal atoms. Since spin-dependent amplitude grows with 
nuclear charge Z faster than Z 2 , one can expect that 
for Fr this amplitude will be significantly larger than for 
other alkali-metal atoms. Besides that, a large number of 
odd isotopes with non-zero nuclear spin makes it possi- 
ble (at least in principle) to study dependence of nuclear 
spin-dependent amplitude on the nuclear structure. 

b. Theory. It is known that parity-nonconserving 
electron-nuclear interaction can be divided into two 
parts: nuclear spin-independent part and nuclear spin- 
dependent one. The respective PNC Hamiltonian can be 
written as follows (atomic units are used throughout the 
paper) : 
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where Gf = 2.2225 x 10~ 14 a.u. is the Fermi constant 
of the weak interaction, Qw is the nuclear weak charge, 
k is the dimensionless coupling constant, a = 707, 7i 
are the Dirac matrices, I is the nuclear spin (/ = | 
for the isotope 211 Fr), and p(r) is the nuclear density 
distribution. 

As we mentioned above, there are three main contri- 
butions to the coupling constant k in the spin-dependent 



part of the PNC interaction 



(_l)/+i 



h-l- 



1 
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where the anapole contribution is given by the constant 
K a f3j (I is the orbital angular momentum of the un- 
paired nucleon) , the constant k,i corresponds to the spin- 
dependent weak neutral currents and the term kq w is 
induced by the interference of the spin-independent PNC 
interaction with the hyperfinc interaction. For heavy nu- 
clei constants K a and kq w are proportional to A 2 / 3 B pi, 
and their ratio depends on the dimensionless constant of 
the weak interaction of the unpaired nucleon with the 
nuclear core g M : 



The root-mean-square charge radius for 211 Fr was mea- 
sured to be r rms = 5.566 fm fll5|| . Using the relation 

r n = a/1 r rms , we find r n = 7.186 fm. 

If \i) and |/) are initial and final atomic states of the 
same nominal parity, then to the lowest nonvanishing or- 
der, the electric dipole transition matrix element (ME) 
is equal to: 



(/K,PNC|J) = Yl 



(f\d q \n) (n\H PNC \i) 



Ei — E n 
(f\H PNC \n)(n\d q \i) 



Ef — En 
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J is the total 
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where /xjv and /i are magnetic moments of the nucleus 
and the valence nucleon, correspondingly. The numerical 
factor q is within the limits 1 < q < 3 (see, e. g., g). For 
the unpaired proton g p ss 7, while for neutron g n sa —2 
H (see also [|J). This estimate shows that for the odd 
isotopes of Fr the anapole contribution dominates in (0) . 
We assume that the nucleus is a uniformly charged 
sphere: 



where \a) = \J a ,F a ,M a ) and F = 
angular momentum. 

In our case the contribution of 77si (see Eq. (Q)) is 
negligible, so we consider only the nuclear spin-dependent 
part of the PNC Hamiltonian. The ME of 77sd can be 
written as follows 



(n\H SD \i) = (-1) /+F < +J < y/l(I + 1)(27 + T)5 FnFi 5 MnMi 

(Jn\\Hsu\\ Ji) , (5) 
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where (J„||F S d||^) = ^j{J n \ho7P(r)\\Ji). 

The ME of the operator d q is given by the following 
expression: 



(f\d q \n) = (-1)*/-"/ ( F l If l q £» ) (-l)^»+^+ 1 v /(2F n + l)(2F / + l)| J p n f ^ ) | (J f \\d\\J n ). (6) 
Applying the Wigner-Eckart theorem to the PNC amplitude: 

(/KpncN) = (-1)^-^ ( _% f l q j^J^/.-F/lldpNclIJi,^), 
and substituting Eqs. (pi) and (0) in Eq. (Eh we get the following expression for the reduced ME of the PNC amplitude: 

(J f ,F f \\d PNC \\Ji,F^ = y/l(I+l)(2I+l)(2F i + l)(2F f + l)Y t [(-l) Jf - Ji { j" j p. } { p f % )) (?) 

(J f \\d\\n,J n )(n,J n \\H S T>\\Ji) ( _^F f - Fl / Jn Jf 1 1 / J„ -h 1 1 {Jf\\Hsn\\n,J n ){n,J n \\d\\J i ) ' 
En-Ei + ( lj \ I I F f j\F z F f I j En -E S 

Note that for the transition between the hyperfine components of the ground state 7s i, one has Ji = Jf = i, 
Ei = Ef = Eis , and Fi = Ff — 1 = 1— \. That leads to some simplification of Eq. (0) : 
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where the sum runs over the states of odd-parity with angular momenta J n — i,§. Novikov and Khriplovich 



pointed out, that for alkali atoms contribution of the in- 
termediate states with J n — | is strongly suppressed. If 
these states are excluded from the sum in Eq. (§h, it can 
be further simplified to the form: 
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Eq . (ft) was used in the semiempirical calculations |1G, 
0|, but here we use the more accurate expression <M). 

c. Method of calculations and results. The Dirac- 
Fock-Breit equations were solved self-consistently on a ra- 
dial grid for the core electrons [Is,.. .,6^3]. Then, the va- 
lence orbitals 7s, 7p, 8s, 8p, 9p were constructed in V^ -1 
approximation. The basis set used in calculations in- 
cluded also virtual orbitals up to 25s, 25p, 24d, 15/, and 
I5g formed with the help of the method described in 
M.g§H. 

To find the nuclear spin-dependent PNC amplitude de- 
fined by Eq. (|8j), one needs to sum over intermediate 
states or solve the corresponding inhomogeneous equa- 
tions (Sternheimer [M or Dalgarno-Lewis |2l|] method). 
Here we apply the Sternheimer-Dalgarno-Lewis method 
to the valence part of the problem as described in |1^, g2| . 
Solving inhomogeneous equation we find the answer in 
the Dirac-Fock approximation: 

(7s, F = 5||dp NC ||7s, F = 4) DF = -0.42x 10" 10 i k a.u.. 

(10) 
It is known that core-valence correlations usually play 
an important role for the PNC amplitudes. We first 
solved the random-phase approximation (RPA) equa- 
tions, summing a certain sequence of many-body dia- 
grams to all orders for both operators in the right hand 
side of Eq. (0). Note that after the RPA equations are 
solved for the operator -ffsD, the MEs (ns||iZsD||nf>3) 
are no longer equal to zero. As a result, the interme- 
diate np3 states also contribute to the spin-dependent 
PNC amplitude. We found that their contribution to 
(7s 7 F = 5||d PNC ||7s,F = 4) is about 10%. That contri- 
bution is neglected in the approximation (Q). The RPA 
correction changes the PNC amplitude to: 

(7s, F = 5||dpNc||7s, ,F = 4) RPA = -0.48xlO~ 10 iK a.u.. 

(11 ) 
The core polarization was taken into account by many- 
body perturbation theory. We completely accounted for 
the second order of perturbation theory and partly for 
the higher orders. In particular, we calculated the struc- 
tural radiation and normalization corrections to the PNC 
amplitude. 

Finally, taking into account that the initial and final 
states are the many-electron states one needs to account 



for the excitations from the npj shells (n=2— 6). Respec- 
tive "core" contribution to the spin-dependent amplitude 
was estimated to be —3.5% 

TABLE I: Nuclear spin-dependent PNC amplitude (7s, F = 
5||dpNc||7s,.F = 4) in the units i x 10" 10 k. The 1-st column 
present the result obtained in the Dirac-Fock approximation 
for the Coulomb-Gaunt potential. Following columns present 
corrections discussed in the text. In the column MBPT the 
Brueckner, structural radiation and normalization corrections 
are summed together. In the column "core" contribution of 
the core excitations is given. 



DFB +RPA 
-0.418 -0.058 



+MBPT 
-0.033 



+ "core" Total 
+0.018 -0.491 



All mentioned corrections are presented in Table Q. 
Summing them up, we finally obtain: 

(7s,F = 5||dpNc||7s,F = 4) = -0.49 x 10~ 10 iKa.u.. 

(12) 
According to Table Q the MBPT corrections to this am- 
plitude are relatively small. Therefore, we estimate the 
accuracy of our result to be about few percent. 

It is interesting to compare this amplitude to similar 
amplitudes in K and Cs. The amplitude (||) strongly 
depends on the nuclear spin /, which is different for all 
alkalis. Therefore, it is convenient to rewrite it in terms 
of the matrix element of the electron operator a = 2s, 
as it was done in Ref. |ll| : 

(7s, F f \\d PNC \\7s, F t ) = iDKiFfMFi). (13) 

In this form the parameter D only weakly depends on the 
nuclear spin /. Combining Eqs. (p"2|) and (|l3|), we get: 



D= 1Q- 12 



-0.07, for 39 < 41 K fill, 
-1.4, for 133 Cs (Eg, (14) 

-11.0, for 211 Fr, 



where we took into account the differences in definition 
of the coupling constant k |2^] . 

One can see that the constant D for Fr is almost an 
order of magnitude larger, than for Cs. According to 
the Refs. |23|, p3, the ratio of the spin-independent PNC 
amplitudes for optical transitions ns — ► (n + l)s for Fr 
(n = 7) and Cs (n = 6) is close to 20. That factor 
also accounts for the 1.6 times difference of the weak 
charges Qw in the PNC operator (Q) for the two nuclei. 
Because the interaction of the electron with the nuclear 
AM gives the main contribution to the spin-dependent 
PNC amplitude, one can expect that the constant k also 
grows as A 2 / 3 [g| . That can account for the extra factor 
~ 1.4 for the amplitude (O) in Fr when compared to 
that in Cs. 
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